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Monte Carlo results of diffusion coefficients 
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We report on Monte Carlo results of diffusion coefficients of lattice star polymers trapped inside 
a fixed network (de Gennes model), ft is found that our data are in agreement with the Helfand- 
Pearson exponential factor a = 0.29. For the pre-exponential power law exponent we find = 2. 
In contrast to existing theoretical predictions, we find that the number of arms / leads to a pre- 
exponential factor of the form exp(— 0.75/). 



I. INTRODUCTION 



Entangled star polymers have the unique property that the viscosity does not increase with a power law of the 
molecular weight, as found for linear polymers, but instead increases exponentially. This experimental result was first 
reported by Kraus and Gruvera and subsequently -confirmed by many othersa. Theoretical explanations have been 
given based on extensions of the reptation modeloQ. An essential assumption made by adapting the reptation model 
to entangled stars is that the dominating mechanism for translational motion is the retraction of ends of the arms 
along their average contour, with the simultaneous projection of unentangled loops into the surrounding matrix. As 
the arms retract the free energy increases due to the loss of configurational entropy. De Gennesu was the first to 
explore this problem using a lattice model defined in the next section. He calculated the number of walks for each 
tube length and concluded that the time required to retract ane arm completely would increase exponentially with 
the arm length AT. In a subsequent approach Doi and Kuzuuu obtained the time dependence of the star relaxation. 
The Doi-Kuzuu theory predicts that the longest relaxation time r m increases exponentially with the molecular weight 
and has a power law prefactor 
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The power law exponent j3 has been predicted by Doi and KuzuuS to j3 — 3, by Pearson and HelfandB to j3 = 3/2, 
and has been estimated by Needs and Edwardsa using Monte Carlo methods to /3 « 1.9. The exponential factor a 



where q is the lattice coordination number. In the present case of the sinxrale cubic lattice one has q = 6 and thus 
anp = 0.29389. The numerical estimate of a — 0.195 by Needs and EdwardsQ is significantly smaller than this. Since 
during the time T m the star polymer diffuses by a tube diameter, which is in the order of unity in our lattice model, 
one can use the relation r m D ~ 1 between the terminal relaxation time r m and the diffusion coefficient D in order to 
obtain 

D = c N-P e~ aN . (3) 

Using Monte Carlo simulations Needs and Edwards! found j3 = 0.59, a — 0.37 and c = 0.033 for / — 3. 

With respect to the conflicting estimates of (3 and a we have attempted to clarify this situation by simulating again 
the lattice model proposed by de Gennes. The advance in computers combined with a highly optimized multispin 
coding technique allows us to simulate over time scales that are several orders of magnitude longer than previous 
simulationsQ. Moreover, we have addressed the dependence of the diffusion coefficient on the number of arms /. 
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II. MODEL AND SIMULATION TECHNIQUES 



We simulate a model for a polymer in a gel that was introduced by de Gennes. In this modelEm, a star polymer 
with / arms and arm length TV is represented by fN + 1 monomers on a 3D cubic lattice, connected by a sequence 
of fN steps on lattice edges. One elementary move in this model is made by randomly selecting one monomer, and 
attempting to move it. If all steps connecting the monomer are located on the same edge, it will randomly move to 
one of the six possible lattice sites (the site it is already in, and five new ones). The rate of any monomer to move to 
any allowed lattice site is once per time unit. A configuration of a three-legged star polymer in the two-dimensional 
version of this model is illustrated in Fig. 1. In this configuration, monomers 1, 2, 3, 6, 7, 9, and 14 can move to three 
other locations, the other monomers are frozen. Note that all our simulations were done in the three-dimensional 
model. 
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FIG. 1. Two-dimensional version of the lattice model for a star polymer with three arms (f=3) and N=7 steps. 



A straightforward way to store the configuration of the polymer would be to store the coordinates of all the 
monomers. However, a more efficient implementation of the algorithm can be obtained by storing the coordinates of 
the central monomer only, plus for each arm the sequence of steps taken from this center. In one elementary move in 
a direct implementation of the dynamics, randomly one of the fN + 1 monomers is chosen and proposed to move. 

• (a) The central monomer can only move if the first steps of all arms are in the same direction; in this case, these 
steps will be replaced by a randomly chosen direction, and the coordinates of the central monomer arc updated. 

• (b) One of the end monomers can always move, and the last step is replaced by a random direction. 

• (c) Any other monomer can only move if the two steps connected to it are opposites. In this case, these steps 
are replaced by a pair of randomly chosen opposite steps. 

To obtain that each of the 6fN + 6 possible moves is tried with unit rate, one elementary move corresponds to a 
time increment of At = l/(6/iV + 6). 

It turns out that the diffusion coefficient decays exponentially with arm length, requiring a very efficient imple- 
mentation of the algorithm to extract reliable numerical results. We achieve that by using multispin coding, enabling 
us to make about seven million elementary moves per second per processor on a SG 200 workstation—In the past, 
we simulated the repton mptLel proposed by RubinsteinEJ with the same approach in an electric fieldliil, and in the 
absence of an electric fieldliatJ, and de Gennes model for linear polymerai-i The processors are 64-bit, allowing 
for 64 concurrent simulations. For each of the 64 simulations, we store the coordinates of central monomer as 
(x Cl y c , z c y ' ■ ■ ■ (x Cl y c , z c Y &3 \ Suppose that we denote the k th bit of long integer w as w k , then step i in arm j 
of simulation k is stored in the three bits {a[if + j] k , b[if + j] k ,c[if + j] k }; if the step is in the positive x-, y-, or 
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z-direction, this triplet is equal to {1,0,0}, {0,1,0} or {0,0,1}, respectively; the negative directions are represented by 
{0,1,1}, {1,0,1}, and {1,1,0}, respectively. Note that opposite steps are each other's binary complement. The core of 
the program, that proposes moves of monomers other than the central and end monomers simultaneously in all 64 
simulations, can then be written (in the programming language C) as: 

f lip= (a [p] Aa [q] ) & (b [p] Ab [q] ) & (c [p] Ac [q] ) ; 
nf= —flip; 

a[p] = (a[p]&nf ) I ( rnda &f lip) ; 
b [p] = (b [p] &nf ) I ( rndb &f lip) ; 
c[p]=(c[p]&nf ) I ( rndc &f lip) ; 
a[q] = (a[q]&nf ) I ( (-rnda) &f lip) ; 
b [q] = (b [q] &nf ) I ( (-rndb) &f lip) ; 
c[q] = (c[q]&nf ) I ( (~rndc) feflip) ; 

Here, p and q = p + / are steps connected to the same monomer, the symbols A, &, — and | denote respectively 
the exclusive-OR, AND, NOT and OR operations, and the triplet (rnda, rndb, rndc) represents a random direction. 
Similar statements can be written down for the central and end monomers. 



III. THE DIFFUSION CONSTANT 



The star polymers were simulated over a long time t to t (around 3 • 10 8 for the longest arm lengths), and the 
coordinates of the central monomers were written to a file at regular times. Afterwards, a quick estimate for the 
diffusion coefficient from 
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with t\ — t to t an d to = 0; with this quick estimate we determined the time interval At after which the mean square 
displacement of the center equals the arm length N, and performed a better estimate of D using the same equation 
with t\ = to + At, averaged over all to. The results are presented in table |. 

TABLE I. The diffusion coefficients Df of star polymers with / arms and arm lengths N. 
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The first step in our analysis of the data is to demonstrate the predicted exponential dependency on N as given by 
(||). In Fig. 2 a semi- log plot of the diffusion coefficients D as a function of N is presented. We found that the data 
of the scaled diffusion coefficient Dexp(0.75/) leads to a fairly good collapse of the data for large N, independent of 
the exponent (3 and the exponential factor a, which is demonstrated in Fig. 2. At large N the data obey the expected 
exponential behavior with an effective exponential factor of a e ff ~ 0.37, very similar as in previous simulationsQ. As 
a guide to the eye, the broken line in Fig. 2 is ~ exp(— 0.377V). 




0.1 E-6 ' 1 1 1 1 
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N 

FIG. 2. Semi-log plot of the Monte Carlo data of the scaled diffusion coefficients De°- 75/ for 3 < / < 6 and 2 < TV < 30. 
The broken line corresponds to ~ e _0,37JV . 

Assuming this exponential law we have analyzed the pre-exponential .power law of (0). This is presented in Fig. 3. 
In agreement with previous simulationsQ and theoretical considerations!] one observes for shorter chains an effective 
power law of N , given by the broken line in Fig. 3. However, at TV > fO there is a significant deviation from the 
initial behavior. The statistical error, as given in Table I, cannot account for this. 

z i i i i i i i rn i i i i i i i rz 

D exp( 0.75 f + 0.37 N ) 




10° N 10 1 10 2 



FIG. 3. Log-log plot of the Monte Carlo data of the scaled diffusion coefficients De * The broken line corresponds 
to ~ iV~ 0,6 . The numbers at each curve denote the number of arms /. 
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Therefore, we have reexamined the data assuming the exponential factor of Helfand and Pearson (0). The corre- 
sponding analysis is presented in Fig. 4. In this case the data indicate a pre-exponential power law of N~ 2 . Therefore, 
in the limit of large N the data of the diffusion coefficient seem to obey the formula 



D 1 = N~ 2 c - a HpN-Q.75f + 1.7 



(5) 



(the subscript "1" indicates the single star behavior trapped inside a fixed network, in contrast to the case of con- 
centrated solutions of star polymers). It is interesting to note that the exponential term can be approximated by 
cxp[-a ffP iV-0.75(/-2)]. 
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FIG. 4. Log-log plot of the Monte Carlo data of the scaled diffusion coefficients De 0,75 - f + a29389JV . The broken line corresponds 
to ~ iV~ 2 . The numbers at each curve denote the number of arms /. 

The normalized data D/Di are presented in Fig. 5. 




FIG. 5. Log-log plot of the Monte Carlo data of the scaled diffusion coefficients D/D\, where Di is given by Eq. (5). 



It is important to note that the /-dependence of D\ appears as a pre-exponential factor. Previous suggestionsi'El 
assumed D ~ exp(— a(f — 2)N), based on the assumption that the diffusion of the branch point requires simultaneous 
retraction of all but two of the star arms. The latter approach cannot be reconciled with our data. However, there is a 
much better agreement between our data and recent estimates based on experimentsEj which gave D ~ exp(— 
These two latter findings, Monte Carlo and experiments, strongly support earlier considerations by RubinsteintZI who 
predicted a much weaker dependence on / as compared to the prediction exp(— afN). Therefore the experimental 
workll3 and our analysis strongly support the view that in order to move the center of the star, the £g~ 2 arms need not 
be retracted simultaneously. Our data cannot be used to distinguish clearly between Rubinstein'sEZI and our formula 
exp(— 0.75/) as given in Eq.(5). 

IV. SUMMARY AND CONCLUSIONS 

We have reported on Monte Carlo results of diffusion coefficients of lattice star polymers trapped inside a fixed 
network. We found that our data arc in agreement with the Hclfand-Pearson exponential factor a = 0.29. However, 
the pre-exponential power law exponent = 2 is different from existing theories and coincides with the exponent 
for linear chains. In contrast to existing theoretical predictions, we found that the number of arms / leads to a 
pre-exponential factor of the form exp(— 0.75/). 

We have not undertaken to estimate the finite size corrections to the asymptotic formula (^) using our Monte-Carlo 
data. This would be probably quite useful comparing our results with more refined theories proposed recentlylij. 
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